Rate equations for Coulomb blockade with ferromagnetic leads 
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We present a density-matrix rate-equation approach to sequential tunneling through a metal 
particle weakly coupled to ferromagnetic leads. The density-matrix description is able to deal with 
correlations between degenerate many-electron states that the standard rate equation formalism in 
terms of occupation probabilities cannot describe. Our formalism is valid for an arbitrary number 
of electrons on the dot, for an arbitrary angle between the polarization directions of the leads, and 
with or without spin-orbit scattering on the metal particle. Interestingly, we find that the density- 
matrix description may be necessary even for metal particles with unpolarized leads if three or more 
single-electron levels contribute to the transport current and electron-electron interactions in the 
metal particle are described by the 'universal interaction Hamiltonian'. 
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I. INTRODUCTION 

Spin-polarized electron tunneling is essential to spin- 
based electronics^ and nanoscale magnetics based on the 
spin-transfer effect i^* 3 . Whereas tunneling through a sin- 
gle tunnel barrier, either between two ferromagnets or 
between a ferromagnet and a normal metal, has been 
studied since the mid 1970s, 4 the study of spin-polarized 
transport through mesoscopic double tunneling junctions 
is more recent. Double mesoscopic junctions are of inter- 
est because of the small capacitance of the central re- 
gion in between the tunneling junctions, which allows 
electrons to be transported one by one?&> Experiments 
have been reported both for normal metal leads with 
a ferromagnetic island between the junctions^ 7 - and for 
a normal island with spin-polarized junctions&Siifl A 
large number of theoretical works has dealt with these 

„„„„„ 11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30 

In this work, we consider the case of a normal metal 
island with ferromagnetic leads. If the temperature is 
much larger than the tunneling rates onto or off the is- 
land, electron tunneling is sequential. In that case, quan- 
tum mechanical correlations between electrons in differ- 
ent states are lost because of thermal smearing, and a 
simple description in terms of rate equations applies. De- 
pending on whether the temperature is small or large 
compared to the level spacing in the island, these rate 
equations describe the probability to find a certain num- 
ber of electrons on the island^i^SiS or the occupation of 
the electronic states in the island 

Whereas the rate-equation approach was applied 
straightforwardly to ferromagnetic leads with collincar 
polarizationsj^Si 2 ^ application to leads with non-collincar 
polarizations requires a formulation in terms of the den- 
sity matrices of degenerate levels, not the occupations 
of statesi£2i2i (If spin degeneracy on the normal metal 
island is lifted, e.g., by a magnetic field, scalar rate equa- 
tions remain valid despite polarized leads.) There are 
two reasons for this additional complication: First, with 
non-collinear polarizations, no common quantization di- 
rection exists, and one cannot avoid a formulation of the 
problem in which electrons tunnel into superpositions of 
states with different spin projections^ Second, coupling 



to the ferromagnetic leads slightly lifts the spin degen- 
eracy and leads to a slow precession of the spin on the 
dotiSS The use of density matrices instead of occupation 
probabilities in the rate equation formalism allows for 
the inclusion of correlations between different quantum 
statesi^S 7 . Since the temperature is much larger than the 
escape rate to the leads, only correlations between states 
with the same energy need to be taken into account. 

Density-matrix rate equations were first used to de- 
scribe transport through a metal particle (or a quantum 
dot or a single molecule) with spin-polarized tunnel con- 
tacts in a recent paper by Konig and Martinek 2 ^ (see 
also Ref. l30|) . These authors used the Keldysh formalism 
to derive the density-matrix rate equations for a dot in 
which only one level contributes to transport. The pur- 
pose of the present work is to formulate a density-matrix 
rate equation for quantum dots in which many electronic 
levels contribute to transport and to simplify the deriva- 
tion of Ref. 30. The extension to many levels is relevant 
for the analysis of experimental data, since the majority 
of experiments feature high bias voltages at which more 
than one electronic level contributes to the current^ 

A remarkable result of our study is that a formula- 
tion in terms of density-matrix rate equations is not 
only needed for spin-polarized leads with non-collincar 
polarization directions, but that it may also be neces- 
sary for unpolarized leads or for spin-polarized leads with 
collinear polarization directions if the metal island has 
a large dimensionless conductance g. These relatively 
large metal particles or quantum dots have degenerate 
or almost-degenerate many-electron levels. Correlations 
between the degenerate states persist during the time 
an electron occupies the quantum dot and need to be 
accounted for using a description in terms of a density 
matrix. The origin of the degeneracy is that in large-g 
metal grains or quantum dots electronic interactions are 
described by the 'universal interaction Hamiltonian'.- 9 
With this interaction, many-electron levels with three or 
more singly occupied single-electron levels are degener- 
ate if their spin is not maximal. For example, there are 
four degenerate states with three singly occupied levels 
and total spin S = 1/2. A rate equation in terms of 
scalar occupation probabilities only 40 will fail to describe 
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FIG. 1: Metal particle attached to two ferromagnetic leads. 



correlations between these degenerate states. A detailed 
description of this case will be given in Sec. Ill Cl below. 



II. MATRIX RATE-EQUATION FORMALISM 

We consider a metal particle — or a quantum dot or 
a single molecule — that is attached to a number of fer- 
romagnetic leads via tunneling contacts with a conduc- 
tance much smaller than the conductance quantum e 2 /h. 
A schematic drawing of a metal particle with two leads 
is shown in Fig. ^ in our formulation of the problem, 
we assume that all leads are fully polarized; a partially 
polarized lead is simply represented by two fully polar- 
ized leads with different densities of states and different 
tunneling rates. We assume that the temperature T is 
much larger than the tunneling rates to and from the 
leads. This is the regime for which rate equations have 
been shown to be a valid description of metal particles 
without spin-polarized leads. 



on the charging energy and exchange interaction of the 
metal particle, the voltages on nearby gates, etc. For 
occupation by zero or two electrons, the many-electron 
level is non-degenerate, and we can use scalars po and p2 
to describe the probability to find the particle in a state 
with zero or two electrons, respectively. If the level is 
occupied by one electron only, one needs to use a 2 x 2 
density matrix p\ to fully describe the state of the parti- 
cle. Conservation of probability implies 

Po +trpx +p 2 = 1. (1) 
Without tunneling to and from the reservoirs, pa, p±, 
and P2 are time independent. The time dependence of po, 
pi, and P2 then arises from tunneling of electrons onto 
or off the metal particle: real tunneling processes shift 
the number of electrons on the metal particle, whereas 
virtual tunneling processes cause a precession of the spin 
if the level is singly occupied. The net tunneling rate to 
and from lead a depends on the direction of the polar- 
ization in that lead, which we describe by means of the 
spinors m a (fh a ) pointing parallel (anti-parallel) to the 
polarization direction of lead a, the tunneling rate T a 
for electrons with spin m a , and the distribution function 
f a in lead a. In order to describe both the virtual and 
real tunneling processes, we combine the rate T a and the 
spinors m a and m a into the spinor tunneling amplitudes 



la 



- pl/2 



IV m a , 7 C 



- Pl/2, 



(2) 



with 7i7 Q = 7l7 Q = r Q and 7^ = 0. 

Virtual tunneling processes can be described by the 
effective Hamiltonian^i 



A. Single doubly degenerate level 

In order to make the connection to previous 
wor i cs? i6 i 29 4 30 we jj rg ^ develop the formalism for the case 

of nonlinear transport through a single level. Our ap- 
proach is closely connected to the works by Nazarov 3 ^ 
and Gurvitz^ who used rate equations to describe high- 
bias transport through a sequence of tunnel barriers. For 
a metal particle in which only one level is relevant for 
transport, we need to consider occupation of the level by 
zero, one, or two electrons, at energies £0, £ij and £2, 
respectively. The precise value of these energies depends 



Hi = ^£(l-2/ a (0) 
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£1 — £0 — £ £2 — £\ — £, 



(3) 



where P denotes the Cauchy principal value. Note that 
if £2 — £\ = £1 — £0, one has Hi proportional to the 
unit matrix in spin space and thus [Hi, pi] — 0: virtual 
excitations do not cause a spin precession without inter- 
actions j2£ The time evolution of the scalars po and P2 and 
the 2x2 matrix pi is described by the equations 29 ' 30 
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dpi 
dt 



a a 

rOi#i - Hipi) + Y /«( £ i - £o)7aPo7a - 2 X^ 1 ~ f a ( £l ~ e °)) [ 

a a 



laliPl + Pl7a7a] 



(4) 



(5) 



3 



a a 

whereas the current through each of the tunnel contacts is calculated as^ 4 ^ 

la = - £0)7aP07a - C 1 - - £ o))7iPl7a + /a(E2 - ElMpiTa - C 1 - /a( £ 2 -£l))7iP27 a - (7) 
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For f a = or f a = 1, Eqs. J3J -0 follow from con- 
sidering the escape of electrons or holes from the metal 
particle into vacuum i^SiSL^ The factors f a and (1 — f a ), 
which also appear in the scalar rate equations? 3 ^*^ are 
inserted to reflect the modification of tunneling rates by 
the electron distribution in the leads. This simple way of 
accounting for the presence of electrons in the leads is no 
longer valid when correlations between the electrons in 
the leads and in the metal particle are formed, such as is 
the case in the Kondo effect, 22 i 23 i 24 i 25 i 26 i 43 i 44 We remark 
that, since Eqs. are meant to describe transport 

at the lowest order in T a only, the energy shift implied 
by the Hamiltonian Hi need not be included in the ar- 
gument of the distribution function f a . For that reason, 
Eqs. I@J-(JI|) describe both the case of low bias and high 
bias, in contrast to the formalism of Refs. which is 

appropriate for high bias only. Also note that Eqs. 
© are consistent with probability conservation, Eq. 
and that they reduce to the standard rate equations 3 ^ 3 ^ 
once the polarizations in the leads are collinear. 

B. General formalism 

For the general description, wc consider a normal- 
metal particle with many-electron levels e k , each of which 
has degeneracy j k ■ The many-electron states are labeled 
|fc,m), where m — 1, . . . , j k . The number of electrons in 
state \k,m) is denoted by N k . Real and virtual transi- 
tions between many-electron states are possible because 
of tunneling of electrons between the metal particle and 
the source and drain reservoirs. As before, we assume 
that this tunneling rate is small in comparison to the 



spacing between electronic levels and temperature. In 
that case, we may describe the state of the dot by a set 
of density matrices p k for each many-electron level, and 
we can neglect correlations between states with different 
energy. 

The tunneling Hamiltonian describing the coupling 
of the metal particle to lead a is determined by the 
jk x jk' matrix v a . k k , containing the matrix elements be- 
tween the many-electron multiplets |fc, •) and \k' , ■) with 
Nk = Nk' ±1. In order to make contact to the rate 
equations derived above, we define the jk X jy matrix of 
transition amplitudes k , = {2nv a /'h) 1 / 2 v : ^. k k ,, where 
v a is the density of states of lead a. We define j^. k k , = 
if N k £ N k , ± 1. One may write 7 ±. fefc , = {T a fl*w±^ k , , 
where w^. k k , is dimensionless and r a is the tunneling 
rate through contact a if the metal island is replaced by 
an electron reservoir. For point contacts, the magnitude 
of w a . k k , is set by the value of a wavefunction at the 
location of the contact ml If the degeneracy of the multi- 
plets |fc, •) and |fc', •) arises from angular momentum, the 
matrix structure of w^. k k , is set by the Clebsch- Gordon 
coefficients. With this notation, virtual excitations lead 
to the effective Hamiltonian H k for the multiplet \k, •), 

H * = ^ P J de£(i-2/ a (0) 

a, A;' 

^ 7a;fc,/c / 7a;fc / ,A: *Y a\k,k'^ a\k' ,k /g\ 

£fc - Eft' - £ Sk> - Sk - £, 

Then the appropriate generalization of the rate equations 
(QJ-© and the current formula (7J is 



dpk 
dt 



^{p k H k - H kPk ) + [f a ( £ k - Efc')7a;fc,fc'P*'7„ i v,* + (! - - £ k)h a -,kM'Pk'^tk'.k 

a, A:' 

- 2 X] [f<*( £ k> _ £ k){la;Kk'lXk',kPk + Pkla;k,k'lt,k',k) 
Ca.k 

+ (1 - f a (£k - £k'))hi, Kk ,l~k',kPk + Pkltk,k>^k> ,*:) 

C Y1 [f°<( £k - £ k'hi;k : k'Pk'la;k',k ~ (! " Z"( £ fc' ~ £ k)ha-M : k'Pk'lXk',k 



k,k' 



(9) 
(10) 



Here the summation over k' extends over all many- electron states different from k. One easily verifies that 
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Eq. conserves the total probability J^k ^ r Pk = 1- 

C. Unpolarized leads 

The density-matrix rate equations @ and (|10|) do not 
necessarily reduce to the standard scalar rate equations 
of Refs. 1.3 III.35I when all leads are unpolarized or when the 
leads have collinear spin polarizations. The reason for 
this, at first, surprising fact is that overlaps between dif- 
ferent many-body states are not accurately described by 
scalar transition probabilities if there are degeneracies. 
With degeneracies, non-orthogonal coherent superposi- 
tions of many-electron states are involved in the trans- 
port process. 

Although level repulsion rules out degeneracies in the 
single-particle states in a generic metal grain or quan- 
tum dot, for large metal grains or quantum dots, (near) 
degeneracies may occur in the many-electron spectrum. 
The origin of the degeneracy is that electron-electron in- 
teractions in metal particles or quantum dots with large 
dimensionless conductance g are described by the 'uni- 
versal interaction Hamiltonian'^ 

H CC = E C N 2 + JS 2 , (11) 

where N is the total number of electrons on the metal 
particle, Eq is the charging energy, S is the total spin, 
and J is the exchange interaction strength. According to 
Eq. (|llfl . the energy of a many-electron state depends on 
the occupation of the single-electron states and the total 
spin S only. This gives rise to degeneracies in many- 
electron states with three or more singly occupied single- 
electron levels. For example, in a metal grain with single- 
electron levels labeled 1, 2, and 3, the two states 




+c} lC { 3 4 3 |0)) , (12b) 

both have three singly occupied single-electron levels 
with total spin 5 = 1/2 and S z = —1/2. Hence, ac- 
cording to the 'universal interaction Hamiltonian', they 
are degenerate. Since both states have the same value 
of S z , the degeneracy is not broken by a magnetic field. 
However, in principle it may be lifted by non-universal 
residual interactions that are not included in the 'univer- 
sal interaction Hamiltonian'^ but such residual interac- 
tions are weak if g ^> 1, and they can be neglected if the 
level splitting that they cause is smaller than the level 
broadening due to escape through the tunnel contacts. 
The degeneracy may also be lifted in metal particles with 
spin-orbit scattering if the spin-orbit rate h/r so is com- 
parable to the mean spacing A between single-electron 
levelsiS^ 



We now illustrate how this degeneracy necessitates the 
use of a matrix rate equation using the example of a metal 
particle with three spin-degenerate single-electron levels. 
For the ease of argument, a magnetic field is applied along 
the negative z axis. We consider transitions from the 
three two-electron states with S = 1, S z = —1, 

|1) =4 2 c{ 3 |0), (13a) 
|2) = |0) , (13b) 

|3> = c\ lC \ 2 |0) , (13c) 

to the degenerate three-electron states (|12|l . As the states 
()13|l are non-degenerate, they are described by means 
of the probability pj of finding the system in state 
j = 1, 2, 3. On the other hand, the states 112|l are degen- 
erate and we need to describe their occupation by a 2 x 2 
density matrix p, 

'-(££)■ <"> 

Transitions from the states (|13J1 to the doublet 1)12(1 occur 
at rates Tj, j = 1, 2, 3. Writing down the time evolution 
of p that results from those transitions, we find 

dp T 1 (I e¥ \ T 2 / 1 e ^ \ 

d7 = y pi U^ 1 ) + Y P2 {e^ 1 ) 

;)+•- as) 

where the remaining terms describe processes that do 
not depend on the pj, j — 1,2,3. Clearly, there is no 
basis that would diagonalize all three matrices in Eq. i|15|) 
simultaneously for arbitrary choice of the pj . This proves 
that it is imperative to use the full matrix structure for 
p in order to properly deal with correlations between the 
states (|T2l . 



III. APPLICATION TO SPIN-POLARIZED 
TRANSPORT 

We now apply the formalism of the previous section 
to transport through metal particles with ferromagnetic 
contacts. We first consider the simpler case of a metal 
particle in which only one energy level participates in 
transport, and then consider the case of multiple levels. 



A. Single doubly occupied level 

The linear-response conductance G of a metal parti- 
cle coupled to two fully polarized ferromagnetic leads, 
labeled L and R, is easily calculated from Eqs. QJ-Q), 
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G = G cos 2 (6>/2) 



1 



4a 2 T L r R sm 2 (9/2) 



[a 2 + (1 - f( £l - £o ) + f(e 2 - e 1 )) 2 ](T L + r R )\ 
where 9 is the angle between the polarizations of the ferromagnets, Go is the linear conductance for 9 — * 0, 

e 2 T L T R (1 - f(e 2 - e 1 ))j(e 1 - e )(l - /(ei - e ) + /fo - £i)) 



Go 



and 



a = P 



(r L +r fl )(H-/( ei -eo)-/(ea-ei)) 
rfC (l-2/(0)(e 2 + eo-2ei) 



2?r (ei - E - £)(e 2 - £\ - £) ' 



(16) 



(17) 



(18) 



For partially polarized leads with polarization P a = (T a — T a )/(r Q + T a ), where T a is the tunneling rate for 
electrons with spin fh a , one finds 



G 



2G 
D 



T L (1-P L )+T R (1-P R ) 



p 2 p 2 r L r R sm 2 > 



T L (l + P R )+T R (l + P L )-DE 



(19) 



where Go is given by Eq. I|17l) above, and we abbreviated 

D = T R (l + P L )(1 - P R ) + r L (l - P L )(l + P R ) + AT L r R P L P R sm 2 (9/2)/(T L + T R ), 

E = a 2 {\ + P L )(1 + P R )(T L + T R )[a 2 + (1 - f{e l - e ) + f(e 2 - £ 1 )) 2 ]- 1 [ri(l + Pr) + T R (l + P L ) 



For symmetric contacts, Yl = Y R and Pl = Pr , 
Eqs. I|16(l and l|19|) were previously obtained in Refs. l29l 
l30l Without the spin-precession term (the first term on 
the r.h.s. of Eq. J5J)), we recover the linear conductance 
calculated by Usaj and Barangerji& after correction of a 
technical mistake in Rcf. 

As pointed out by Konig and Martinek, the role of 
the spin-precession term is to reduce the angular depen- 
dence of the conductance. Our general results (|16fl and 
(I19J1 show how this reduction depends on the symmetry 
of the contacts: the reduction is strongest for symmetric 
contacts (T R = T^), whereas it vanishes in the generic 
case of very asymmetric contacts (r^ < Ti). In the 
latter case, the spin precession axis is aligned with mi; 
precession around rriL does not change the angular de- 
pendence of the conductance. 



B. Case of up to three electrons on the dot 

For a metal particle in which more than one level con- 
tributes to transport we calculated the differential con- 
ductance G = dl/dV numerically as a function of the 
bias voltage V. 

The numerical calculation was done for a metal par- 
ticle in which five single-electron levels, with a total of 
two or three electrons, contributed to the current. The 
lead polarizations were chosen parallel or anti-parallel, 
with polarizations Pl=Pr = P- With a maximum of 
three singly-occupied levels, the largest possible spin on 
the dot was 3/2. The positions of the single-electron 



energy levels were taken from the center of a matrix 
drawn from the Gaussian Orthogonal Ensemble of ran- 
dom matrix theory, and the temperature T was set at 
one percent of the mean spacing A between the single- 
electron levels, to ensure that all features in the current- 
voltage characteristic could be resolved in the numeri- 
cal calculation. Electron-electron interactions were de- 
scribed using Eq. (JTTJ. In the numerical calculations, 
we set E c = 25A, and J = -0.32A. (Values of the 
exchange constant J are tabulated in Ref. f° r most 
normal metals.) The tunneling rates Yl and T R were 
chosen < 0.1k R T and equal for all levels, as is appropri- 
ate for metal particles with wide tunnel barriers. The 
source-drain voltage V was applied to the right lead and 
was assumed to change the effective chemical potential 
in the right lead only. 

The use of leads with collinear polarizations in the nu- 
merical calculations eliminates most of the necessity of 
using density-matrix rate equations, except for the de- 
generate S = 1/2 states with three singly-occupied levels. 
These states are fourfold degenerate. We denote them 

1^. +> = ITU) + e"^ |UT) + |m», (20a) 

\\,-) = -^( e=¥± ITU) + ^ IT-IT) + |iTT», (20b) 

I - = -±=(e^ lUD+e^ |m) + |UT», (21a) 

\-\,-) = 4|( e ^ ITU) + ^ I4U) + |UT»- (2ib) 
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However, only the twofold degeneracy inside the pairs 
with S z = 1/2 and S z = —1/2 is relevant, and it is 
sufficient to describe the occupation of the four 5 = 1/2 
states with two 2x2 density matrices p(S z = 1/2) and 
p(S z = -1/2) as in Eq. (JUJ. 

To illustrate the use of the matrix rate equations in 
this case, we write down the full transition vectors for 
the transition between the S = 1 triplet states and the 
two doublets and (2TJ. We denote the S = 1 triplet 
state by \S Z ), with S z = —1,0,1. In relation to the 
two energy levels already occupied in the triplet state, 
we consider adding an electron in a single-electron level 
with higher, lower, or intermediate energy and denote 
the different vectors by superscripts h, I, and m, respec- 
tively. Choosing the orientation of the leads as the spin 
quantization axis, the nonzero transition vectors for ad- 
dition of an electron with spin up from lead a as they 
appear in the rate equation (JSJ are Ja' h 



7, 



y = (r» a )Va«,+.l, and 
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FIG. 2: Probability of finding total spin 3/2 on the dot for 
anti-parallel (top panel) and parallel lead polarizations (bot- 
tom panel) featuring polarizations P—Pl~Pr=0.95, 0.85, 
0.75, 0.65, 0.55 (top to bottom), with F R /F L = 0.2. 



The transition vectors for adding a spin-down electron 
follow straightforwardly from the above. The amplitudes 
for removing electrons are obtained from the above by 
hermitian conjugation. Although these were not consid- 
ered in the numerical calculations, we mention that the 
overlap matrices for non-collinear lead polarizations can 
be obtained from Eqs. (122H by combining the transition 
vectors into 4x3 matrix amplitudes for the full transi- 
tion from the spin-1 triplet to the three-electron spin-i 
quadruplet, followed by multiplication with appropriate 
representations of rotation matrices. In this particular 
case, the transition matrix amplitude would have to be 
multiplied with a four-dimensional representation from 
the left and a three-dimensional representation from the 
right. The relevant rotation matrices are listed in the 
appendix. 

We now turn to the results of our numerical calcula- 
tions. One expects that anti-parallel lead polarizations 
cause spin accumulation on the metal island. This can in- 
deed be observed in our solution of the rate equations, as 
shown in Fig. [21 where we plot the probability of finding 
spin 5 = 3/2 (as opposed to S = 1/2) on the metal par- 
ticle for different lead polarizations P. For anti-parallel 



lead polarizations, the probability to find S = 3/2 in- 
creases with increasing polarization, whereas it is virtu- 
ally independent of polarization for parallel lead polar- 
izations. 

In Figs. |3 and 0] we address the dependence of peaks 
of the differential conductance on the lead polarization 
P. The figure describes both parallel lead polarizations 
(positive P in the figure) and anti-parallel lead polariza- 
tions (negative P in the figure). Some peaks evolve non- 
monotonously, whereas others rise or fall monotonously 
in magnitude. Of particular interest are those conduc- 
tance peaks that evolve from positive to negative values. 
Negative conductance peaks at voltage V arise if, upon 
reaching that voltage, a many-body state that is poorly 
connected to other states is made accessible. (This type 
of behavior is not limited to ferromagnetic leads.) Not 
only does there seem to be a tendency toward negative 
differential conductance upon going from anti-parallel to 
parallel leads but also when making the transition from 
<C to a more symmetric coupling ~ Tl, as 
shown in Fig. 01 In general, conductance peak spectra 
for a different polarization orientation can therefore look 
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FIG. 3: Excerpt from the spectrum of several conductance 
peaks for anti-parallel and parallel lead orientations as a func- 
tion of source-drain voltage V (in units of the mean level 
spacing A) and of lead polarization P = Pl — Pr. For ease 
of presentation, the case of parallel polarization is plotted 
against negative polarization. The excerpt shown here does 
not include the dominant low-energy peak. 
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FIG. 4: Normalized differential conductance as a function 
of source-drain voltage V for parallel (left panels) and anti- 
parallel (right panels) lead polarizations P=Pl=Pr=0.9. 
The vertical scale has been normalized to the magnitude G max 
of the largest (first) conductance peak. 



very different in terms of both position and magnitude of 
the most dominant peaks, even if recorded for the same 
sample. 

IV. INFLUENCE OF SPIN-ORBIT 
SCATTERING 

The role of spin-orbit scattering inside the metal parti- 
cle is best illustrated by considering the case of transport 
through one doubly (Kramers) degenerate level. (Spin- 
orbit scattering in the reservoirs instead of in the metal 
island was considered in Ref. 00 With spin-orbit scat- 
tering, the two eigenfunctions of the level are spinor wave- 
functions. Once the spin quantization axes are fixed at 
a reference point in the metal particle, the two spinor 
wavefunctions define a spatially dependent spin quanti- 
zation axis in the metal particle. The spinors m a and rh a 
that define directions parallel and anti-parallel to the po- 
larization direction in lead a are defined with respect to 
the quantization axis at the contact with lead a. Hence, 
the presence of spin orbit scattering in the normal metal 
particle alters the spinor structure of the transition am- 
plitudes 7c , but it does not change the general structure 
of the rate equations. The same conclusions hold for the 
general case. 

The above considerations imply that, for a metal par- 
ticle coupled to ferromagnetic source and drain reservoirs 
via two point contacts and with only a single level con- 
tributing to transport, the sole effect of spin-orbit scat- 
tering is a sample-specific shift of the angle between the 
polarizations in the two leads. On the other hand, for 
a metal particle coupled to source and drain leads via 
many-channel tunneling contacts, or for a metal particle 
in which more than one level contributes to transport, the 
effect of spin-orbit scattering is more complicated since 
different levels and different channels experience different 
rotations of the spin reference frame. In particular, with 
strong spin-orbit scattering, (spin) transport through dif- 
ferent channels or through different levels will involve 
completely different rotation angles, so that the effec- 
tive degree of spin polarization in the junction is greatly 
reduced. In the limit of a large number of channels and 
strong spin-orbit scattering, the rate equations in fact 
reduce to the unpolarized case. 



V. CONCLUSION 

We have extended the rate-equation formalism to the 
case of a normal island (metal nanoparticle, quantum 
dot, or single molecule) attached to spin-polarized con- 
tacts with non-collinear polarization directions. Our for- 
malism provides a transparent description of the sequen- 
tial tunneling process in this system, and is suitable for 
applications to both linear and nonlinear transport. 

Whether one has to employ matrix rate equations or 
the simpler scalar rate equations is determined by the 
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symmetries and energy degeneracies of the metal island 
and the leads. We distinguish the following cases: (i) 
leads and island have the same symmetries and degen- 
eracies, (ii) some symmetries present on the island are 
broken in the leads, and (iii) leads and island feature 
the same symmetries but there are additional degenera- 
cies on the dot. Case (i) corresponds, e.g., to a normal- 
metal island with unpolarized leads and spin-degenerate 
energy levels on the dot, or to a normal- metal metal 
island with spin-polarized leads with collincar polariza- 
tion directions. In that case, the standard rate equations 
of Refs. I34ll35l are applicable if there are no further de- 
generacies on the island. Case (i) also describes spin- 
polarized leads with non-collincar polarization directions 
if spin degeneracy on the island is lifted by a magnetic 
field, be it an applied field or the stray field of the fer- 
romagnetic leads. In contrast, the other two situations 
require matrix rate equations: Case (ii) applies to normal 
metal particles with ferromagnetic leads that are polar- 
ized in non-collinear directions so that tunneling occurs 
into coherent superpositions of states on the dot. The 
additional degeneracies required for case (iii) arise in the 
many-electron spectrum of generic metal particles with or 
without spin-polarized leads if the 'universal interaction 
Hamiltonian' describes the electron-electron interactions 
on the metal particle. Such degeneracies can be lifted by 
non-universal interaction corrections if the resulting en- 
ergy splitting is larger than the level broadening due to 
escape to the leads. These corrections scale as l/g, where 
g is the dimensionless conductance of the metal particle, 
and as a consequence, the larger the size of the normal- 
metal island, the more important the use of matrix rate 
equations becomes. 
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APPENDIX: ROTATION MATRICES 

This appendix contains the rotation matrices for trans- 
formation from the basis aligned with the quantization 
axis to a basis forming a relative angle 9 with the orig- 
inal quantization axis. We restrict ourselves to the case 
of up to three electrons on the dot, so that we only need 
representations up to dimensionality four. (The four- 
dimensional representations correspond to the quadru- 
plets of total spin S = 3/2 and S = 1/2 in the case 
of three singly-occupied single-electron levels.) Also, we 
only consider polarizations in the xz plane, so that any 
superposition of single spins can be expressed in terms of 
real coefficients. 

1. Two dimensions: spin-1/2 doublet 

In the standard basis of spin up and spin down, 

(1,0) = It), (o,i) = U), (A.i) 

the rotation matrix is 

R 2 = ( C ° S{9 ( 1 S[n %\ ) . (A.2) 
V — sm(fc') cos(#) J y ' 

1. Three dimensions: spin-1 triplet 

In the basis created by successively applying the spin- 
lowering operator S- to 

(1,0,0) ee |TT) , 

(0,1,0) ee (|T|) + ||T»/\/2, (A.3) 

(0,0,1) EE |||>, 

the rotation matrix reads 



cos 2 (6>) v / 2cos(6»)sin(6») _ sin 2 (6>) 

Ra=\ -V2cos(6>)sin((9) cos 2 (<9) - sm 2 (0) ^2 cos(0) sm(9) |. (A.4) 
sin 2 (<9) -y2cos(6»)sin(6») cos 2 (6») 

3. Four dimensions: spin-3/2 quadruplet 

The maximum spin-3/2 states feature electrons in three singly occupied single-electron levels. In the basis obtained 
by successively applying the spin-lowering operator S- to |TTT), 

(1,0,0,0) ee |TTT) , (0,1,0,0) ee (mi) + |UT) + liTT»/V3, 

(o,o,i,o) = (|TU> + im> + IUt))/V3, (0,0, 0,1) ee mi), (A.5) 
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we obtain the following rotation matrix 



R 



4-t 



/ _cos 3 (0) V3cos 2 (6»)sin(6») V3 cos(0) sin 2 (<9) sin 3 (6») ^ 

-\/3 cos 2 (<9) sin (<9) cos 3 (<9) - 2 cos(6>) sin 2 (<9) - sin 3 (<9) + 2cos 2 (6>) sin(0) v / 3cos(6») sin 2 (0) 

\/3cos((9)sin 2 (60 sin 3 (6>j )- 2 cos 2 (0) sin(0) cos 3 (6>) - 2 cos(6>) sin 2 (<9) V3 cos 2 (6») sin(0) 

\ -sin 3 (0) Vtcos(0)sin 2 (6») -V3 cos 2 (0) sin(0) cos 3 (0) J 



(A.6) 



4. Four dimensions: spin-1/2 quadruplet 

For the four- fold degenerate many-electron state with total spin 5 = 1/2 and three singly occupied single-electron 
levels, we write the rotation matrix in the basis 

(1,0,0,0) = (e 2 ^ ITU) HIT) + |ITT))/V3, (0,1,0,0) = [e=^ \\\\) + |UT) + IITT))/V3, 

(0, 0, 1, 0) = (e 2 ^ ITU) + (IT!) + |UT))/\/3, (0, 0, 0, 1) = (e=^ \Ui) + ^ \\\\) + |UT»/>/3. (A.7) 

The first two vectors have = 1/2, the other two have S z = —1/2. In this basis, the rotation matrix then reads 



2 



cos (9) 

cos (0) 
| sin (0) e 2 -^ 

i, 



• sin (9) e $ 



■\sm{6)e^ i \ 

-\sm{6)e = ¥ L 

' cos(0) 

cos (9) 



(A.8) 
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